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Abstract 

A condensate of charged scalars in a neutralizing background of fermions 
(e.g., condensed lielium-4 nuclei in an electron background in white dwarf 
cores) is investigated further. We discuss an effective Lagrangian approach to 
this system and show that the strong screening of an electric charge found pre- 
viously in arXiv:0806.3692 in a mean-field approximation, is a consequence of 
a cancellation due to a phonon. The resulting propagators contain terms that 
strongly modify their infrared behavior. Furthermore, we evaluate a one-loop 
fermion quantum correction to the screened potential, and find that it is also 
suppressed by the phonon subtraction. Therefore, charged impurities (e.g., 
hydrogen or helium-3 nuclei) will be screened efficiently by the condensate. 



1. Introduction and Summary 



Consider a system of negatively charged fermions (e.g., electrons) of mass mf and 
of number-density Jq, and positively charged scalars (e.g., helium-4 nuclei) of mass 
niH, with the net charge of the system being zero. At densities Jo ~ (0.5 — 5 MeV)^ 
the average inter-fermion separation, li ~ Jq ~ (400 — 40 /m), is much smaller 
than the atomic scale ~ 10^ /m, while being much greater than the nuclear scale 
~ 1 fm. Therefore, neither atomic nor nuclear effects would matter at such densities. 
The fermions would form a degenerate gas with the Fermi energy exceeding inter- 
fermion repulsion energy, which can be ignored. These are the conditions thought 
to be present in cores of helium white dwarfs. 

As a helium dwarf star cooled from ~ 10^ — 10^ down to lower temperatures, 
the hclium-4 nuclei could have formed a crystalline structure. However, the zero- 
point oscillations of helium-4 ions at these densities would "melt" such a crystal, 
and it was proposed [1, 2], that the helium-4 ions could instead condense into a 
macroscopic state of large occupation number - the charged condensate - specific 
properties of which were discussed in detail in [2]. 

In particular, a static potential between probe charges placed in the condensate 
were calculated in [2] in the mean-field approximation, and were shown to be sup- 
pressed by the exponential factor, ~ exp(— Mr)cos(Mr)/r, where the scale M is 
defined in terms of the electric charge e, mass m^^, and number-density Jq: 

M^(2eWo)'/' >JI'^. (1) 

Although this strong screening may well be a reason why the condensation of charged 
bosons takes place in the first place, the following issues emerge in this regard: 

(i) It may seem that the exponent exp(— Mr) is due to a state of mass M. The 
distance scale 1/M is shorter than the average inter-particle separation - an effective 
short- distance cutoff of the low-energy theory. Then, a state of mass M, if existed, 
would have been beyond the scope of the low-energy field theory description, and 
the above-cited potential would have been unreliable. 

(ii) The momentum-space static potential, the Fourier transform of which is 
proportional to ~ exp(— Mr)cos(Mr)/r, reads as follows (see Section 3) 

G{u = 0, k) ^ i^' + + j , (2) 

where, stands for a certain quantity such that \rn^\ ^ (sec Ref. [2] and 
Sections 3,4 below). The term M'^/k'^ gives rise to a significant modification of the 
Green's function in the infrared. The interpretation of such an infrared-sensitive 
term in the denominator of a propagator calls for an explanation. 

(iii) In the niH — ^ oo limit one would expect the heavy scalars (helium-4 nuclei) 
to decouple. It is not exactly clear from (2) how such a decoupling takes place, and 
what is its interpretation. 
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The above three issues, can be clarified, as we will summarize below, and show 
in this paper. The points (i) and (ii) get resolved by observing that the Green's 
function (2) can formally be decomposed as follows: 



The first term on the r.h.s. can be thought of as an instantaneous screened Coulomb 
potential, while the second one can be interpreted as a potential due to a phonon. 
The phonon in this case is a collective excitation of motion of charged scalars within 
the fermion background. As was shown in [2], the phonon is a light excitation 
that belongs to the spectrum of a low-energy effective field theory. It's just the 
cancellation due to this light mode that gives rise to the exponential exp(— Mr), 
and not a hypothetical state of mass M. 

Before we turn to the point (iii) let us make two important comments. First, 
wc should note that the second term in (3) has three poles. The residue of the pole 
at + rriQ = cancels exactly with that from the first term of (3). The remaining 
two poles, describe both the heavy state of mass ~ rriH, and a light state, which 
actually is the phonon [2]. For simplicity of the discussions, in this work we'll be 
using a somewhat arbitrary language by calling the whole second term in (3) the 
phonon contribution. 

Second, wc note that the form of (3) may be puzzling itself: the negative sign in 
front of the second term on the r.h.s. may seem to be suggesting that the phonon is 
a ghost. This indeed would have been a ghost contribution if we were to obtain such 
a term in a spectrum of a Lorentz-invariant theory. Here, the background solution 
fixes the Lorentz frame, and the spectrum of small perturbations is not Lorentz- 
invariant [1]. Because of this, the negative sign in front of the phonon Green's 
functions is not a signature of a ghost. In particular, as we'll show in Section 3, 
the Hamiltonian for the fluctuations that give rise to (3) is positive semi-definite! 
Moreover, the measure of Lorentz-breaking is the scale M. In the M — > limit, the 
phonon term in (3) vanishes, as one would expect. 

Let us now turn to the point (iii). In the limit rriH ~^ oo, which implies M — > oo, 
the phonon effects should go away. This is reflected in (3) as vanishing of the 
whole potential, which is just a consequence of using the static approximation, 
and can be understood in the following way: The phonon mixes with the timelike 
component of the gauge fleld, and due to this acquires an instantaneous part. Then, 
the instantaneous parts in (3) cancel between the gauge and photon contributions. 
However, the dynamical part of the phonon is also reducing to zero, because the 
group velocity of the phonon vanishes in the uih oo limit, as we will show in 
Section 3. 

Having the above properties clarifled, in this work we will study further the 
static charge screening in charged condensate. The physical question is that of 

static interactions between charged impurities placed into the condensate. As we 
mentioned above, the scalar condensate gives rise to the suppression of the static 
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(3) 



2 + k2(k2 + m2)2/4M4 ■ 
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potential by a factor exp(— Mr). Inclusion of the fermion fluctuations via the mean- 
field Thomas- Fermi method, does not change the above result significantly [2]. 

In Section 4 we go beyond the mean-field approximation and include a one-loop 
quantum correction due to the polarization diagram. This diagram is suppressed by 
an additional power of the electromagnetic coupling constant ctem = e^/47r, and one 
would expect the quantum correction to be insignificant. 

However, this is not the case for the following subtle reason. The one-loop correc- 
tion introduces branch cuts in the Green's function (3), which give rise to additional 
contributions to the static potential in the position space. These additional terms 
have oscillatory nature with a power-like decaying envelope. Even though they are 
formally suppressed by Ola"^^), to a good approximation they end up being 0{1), 
and can dominate over the exponentially suppressed term at sufficiently large dis- 
tances. 

Such an oscillatory potential is know in non-relativistic fermion systems as the 
Priedel potential (for original references see [3]), which was generalized to relativistic 
fermions by Sivak [4] and by Kapusta and Toimela [5]. Here we perform analogous 
calculations for the charged condensate and find that the phonon subtraction sup- 
presses further the Priedel-like potential by the factor Jq^'^/M^ ~ 10~^^. In spite of 
this, for reasonable separations between the test particles, the Priedel-like potential 
dominates over the exponential one. 

Therefore, charged impurities, such as hydrogen or helium-3 nuclei, once placed 
within the charged condensate, will interact very weakly with the ambient con- 
densate and with each other. The impurities are screened very efficiently by the 
condensate! 

Let us make a few comments on the literature. Condensation of non-relativistic 

charged scalars has a long history, the original works being those by Schafroth [6], 
in the context of superconductivity, and by Foldy [7], in a more general setup. 
An almost-ideal Bose gas approximations was assumed in those studies. For this 
assumption to be true, densities had to be taken high-enough to make the average 
inter-particle separation shorter than the boson Bohr radius [7]. In terms of our 

parameters this would be the case if Jq^^ ^ctcmfT^H- 

However, for the helium-electron system the above condition would translate 
into super-high densities, at which nuclear interactions become significant. Instead, 

1 /3 

in this work we're studying charged condensation in the opposite regime, Jq <C 
o^em'iTT'H, where the nuclear forces play no role. Moreover, our mechanism of screening 
the electric charge is different, as it is based on cancellation by the phonon. 

A possibility of having a changed condensate in helium white dwarfs was previ- 
ously discussed in Ref. [8]^, where the condensation was studied using an approx- 
imate variational quantum-mechanical calculation in conjunction with numerical 

^We thank Andrei Gruzinov for bringing this paper to our attention. Regretfully, this happened 
after our work on charge condensation in white dwarfs [2] was already published. We use this 
opportunity here to emphasize differences between our approach and that of [8]. 
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insights in a strongly-coupled regime of electromagnetic interactions. A degree of 
reliability of such a scheme is hard to assess^. Finally, in the context of a relativistic 
field theory condensation of scalars was discussed in, e.g., Refs. [9, 10, 11], some of 
the results of which we cite below. 

A. Dolgov, A. Lepidi, and G. Piccinelli [12] (appeared on archive the same day 
with this paper) calculated a one-loop correction to the photon sclf-cncrgy in the 
background of charged condensate at nonzero temperature. These author find sim- 
ilar infrared modifications of the static potential as we did. Moreover, their results 
also include nonzero temperature effects, and in that part, are more general than 
ours. On the other hand, our emphasis is on understanding of these results in terms 
of the phonon and effective Lagrangian language. 

2. Effective Lagrangian and its Solution 

In this Section we discuss an effective Lagrangian description of charged conden- 
sation. We will treat the bosonic and fermionic degrees of freedom separately. This 
has two reasons: (a) There is a significant disparity of physical scales in these two 
sectors; (b) Due to the Fermi surface, the dynamics for fermionic fluctuations should 
be though of in terms of near-the-Fermi-surface degrees of freedom. For instance, 
when the energy of these excitations scales to zero, their momenta scale to the Fermi 
surface [13]. To be able to describe the dynamics of the near-the-Fermi-surface fluc- 
tuations, including their one-loop quantum polarization effects, we will adopt for 
them a microscopic Lagrangian with a relativistic chemical potential 

i;{i^''Di,-mf)ip + Hfip+il!, (4) 

where the covariant derivative for fermions (electrons) is defined as Z?^ = dfj^ + ieAfj^, 
and fif denotes the fermion chemical potential, that equals to the Fermi energy at 
zero temperature Hf = ep = [(Stt^Jq)^/^ + m^]^/^ . 

For the bosonic sector, on the other hand, we will employ the effective order- 
parameter description. In particular we will look for classical solutions of the equa- 
tions of motion of the effective order-parameter Lagrangian. 

How could a classical solution describe the condensation which is an inherently 
quantum phenomenon? Denote the particle creation and annihilation operators by 
Qq and ao respectively; then, the quantum-mechanical noncomutativity of these 
operators, Oq oo — oooj ~ K becomes an insignificant effect of order 0{h/N), when 
the number of particles in the condensate state, (o-o '^o) ~ is large enough, N ^ 1. 
Thus, the classical description of the coherent state with a large occupation number 

^Furthermore, using the ordinary neutral Bose-Einstcin (BE) condensation to describe the 
charged condensate, as it is done in a number of works in the hteraturc, is unjustified. Properties 
of a neutral BE condensate differ significantly from those of the charged condensate considered 
here. For instance, specific heat at moderate temperatures in the former is due to a phonon gas, 
while in the latter it is due to the degenerate electrons (see, e.g., [2]). 
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- the condensate - should be vahd to a good accuracy [14]. On the other hand, 
collective excitations of the condensate itself should be quantized in a conventional 
manner. 

The above arguments lead to the following decomposition of the order-parameter 
operator describing the condensate: 

$ = + 5$ , (5) 

where $c« denotes just a classical solution of the corresponding equations of motion, 
and describes the condensate of many zero-momentum particles, while 5$ should 
describe their collective fluctuations. 

In the non-relativistic approximation in mn-i the effective order-parameter La- 
grangian that is consistent with the translational, rotational, Galilean and the global 
U{1) (scalar number) symmetries, as well as with local gauge invariance, can be 
written as follows: 

>Ce// = V (^-{^*Do^ - (Do$)*$) - , (6) 

where Z^o = (<9o ~^2eAo), Dj = {dj — i2eAj), while V{x) stands for a general polyno- 
mial function of its argument. The coefficients of this polynomial are dimensionful 
numbers that are inversely proportional to powers of a short-distance cutoff of the 
effective field theory, V{x) = ''-^"(^"Z^^")- 

Such a Lagrangian was first proposed by Greiter, Wilczek and Witten (GWW) 
[15] in a context of superconductivity. The requirements leading to (6) are the 
following: (a) In the lowest order in fields it gives rise to the standard Schrodinger 
equation for the order parameter^; (b) It respects all the appropriate symmetries of 
the physical system at hand, and gives rise to an operator relation between a current 
of charge q and the momentum density Jj = {q/mH)Toj. (c) It gives an appropriate 
spectrum of Nambu-Goldstone bosons in the decoupling limit^. 

The GWW effective Lagrangian can also describe charged condensation, as we're 
about to show. For this we use the representation of the order parameter in terms 
of its modulus and the phase, (which is well-defined everywhere since the modulus 
will have a nonzero VEV): 

$ = Eexp(ir). (7) 

^In our nomenclature electrons are "particles" while lielium-4 nuclei are "antiparticles" , this 
explains the choice of the plus sing in front of the first term in (6). 

''One could add to the argument of the function P{x) in (6) terms ^nr^*^, A($*<1?)^, etc. The 
former could be absorbed into a gauge potential by a constant shift, but in any case, in the 
condensation phase, which is a point of the primary interest here, the non-relativistic chemical 
potential /zjvh, should be zero (or equivalently, the relativistic chemical potential jjg = (jlnr + fn,H 
should exactly equal to the particle mass /ig = mu)- Furthermore, the existence of the quartic 
term for our system would not play an important role as long as A ;^ 1 and Jo <C m%. 
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Using this decomposition wc rewrite the argument of the function V and the effective 
Lagrangian (6) in the following form: 



(V,-E)2 + (2e)2E2S| 



(8) 



2mH 



where we introduced a gauge-invariant field B^j^ = A^ — d^V/2e. Varying (8) w.r.t. E 
and Aq, and restricting ourselves only to constant fields, we obtain the corresponding 
equations: 



where in the last equation we took into account that the gauge field also couples 
to the background fermion charge density eJo as follows: —cAqJq. Note, that for 
space-time constant fields the equation of motion obtained by variation of (8) w.r.t. 
r is satisfied identically. 

There exists a solution to these equations for which 



Since on the solution the argument of (8) is zero, the condition 'P'(O) = 1 is satisfied 
by any polynomial functions 'P{x) for which the first coefficient is normalized to one 



The above solution describes a neutral system of negatively charged electrons of 
charge density — eJo, and positively charged scalar (helium-4 nuclei) condensate of 
charge density 2e$+$ = 2eE^ = eJo. 

Let us now turn to the issue of fiuctuations about the classical solution. In 
Ref. [15] the Lagrangian (6) was proposed for the description of low-energy theory 
of superconductivity, where a phonon has a linear dispersion relation uj oc |k|, for 
small |k|. Such a relation was obtained in [15] by fixing E(x) to a constant, but 
allowing for fluctuations of the phase V. 

The system considered here is different, however. In the limit of switched-off 
Coulomb interactions (decoupling limit) it should reduce to a collection of free par- 
ticles of mass niH- The latter undergoes Bose-Einstein condensation at low temper- 
atures. The resulting condensate exhibits no super-conductivity/fiuidity because of 
the lack of interactions between particles. Hence, the dispersion relation for the first 
possible excitation in this system should be ~ k^/2mH. 

Remarkably, this dispersion relation can be obtained from the GWW Lagrangian 
(6) by allowing for the S(a;) field to have small fiuctuations. The latter mix with 
the fiuctuations of the phase, and, as a result, gives rise to the dispersion relation 
UJ ~ \s? l2mH, as we will show in the next Section. 




(9) 



2E2 = 



Jo, S^ = 0, P'(0) = 1. 



(10) 



V{x) ^ X + + ... . 



(11) 
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The above arguments suggest that taking a non-relativistic hmit to decouple a 
dynamics of a heavy mode in the presence of a condensate, is a bit subtle, as is 
well known in solid state theory: For instance, even though the ions in a simple 
crystal lattice could be very heavy and non-relativistic, there still could exists a 
low-energy mode, a phonon, with a linear dispersion relation cj ~ |k|, associated 
with a collective motion of those heavy ions. The decoupling of heavy ions in the 
non-relativistic limit manifest itself as vanishing of the group velocity of the phonon 
in the infinite-ion-mass limit. 

A similar effect for a many-body system of free bosons can be sketched as follows: 
Consider a single massive state with the dispersion (^)^ = -|- m\c^. To describe 
a many-body system of such particles one could replace 00 hy uj + /i^c, where /i^ is 
the chemical potential, and put iis equals to mj^, in order to describe condensation. 
Then, the dispersion relation becomes (^)^ = — 2ujmH- Prom this, one gets 
a dispersion relation in the non-relativistic limit: uj — \s? /2mH-i with the group 
velocity Vgr = |k|/m// being a C(l/mj/) effect, consistently with ones expectation 
and the heavy mass decoupling. 

For convenience, in what follows we will start our discussions with a rclativistic 
Lorentz-invariant Lagrangian that retains the heavy scalar mode. Normally, this 
mode would not be included in the effective Lagrangian, and we should be careful to 
separate heavy and light modes in our calculations, and use only the light ones for the 
discussion of physical effects. Thus, the full Lagrangian containing the relativistic 
order-parameter 0, gauge boson, and fermions read as follows: 



A classical nonzero vacuum expectation value of the field 4> can serve as an order 
parameter for the condensation of the helium-4 nuclei, thus describing a state with 
a large occupation number. Fluctuations of the order parameter are expected to 
describe the collective modes of the condensate^. 

The covariant derivative for the scalar reads D^ = d^ — 2ieA^ = — ii^eA^ + 
/'.sf^o) ) where is the relativistic chemical potential for the scalars. The Lagrangian 
(12) is invariant under global f/s(l) transformations, responsible for the conservation 
of the number of scalars. Another global Uf{l) guarantees the fermion number 
conservation. One linear combination of these two symmetries is gauged, and the 
corresponding conserved current is coupled to the photon field^. 

As before, we introduce the representation for the scalar field, = 73"^^'"' 
work in the unitary gauge where the phase of the scalar is set to zero, a — 0. In 

^The Lagrangian (12) is more restrictive than it's needed - the effective Lagrangian (6) allows for 
a more general velocities of propagation than the Lorentz-invariant one. However, this degeneracy 
of (12) will be lifted by the fermion fluctuations, see Sections 3 and 4. 

^We could add the quartic scalar self-interaction term to (12), but this won't change our results 

1 /3 

significantly, as long as the quartic coupling is not strong, and niH ^ Jq ■ 
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this gauge, the Lagrangian density reads: 
>C = -\fI + \{d,af + ^^Ay - ]^mW + i'iil'D, - mj)ij + /./^+^ . (13) 

Since the chemical potential is nonzero, Aq has an expectation value, and this term 
plays the role of a tachyonic mass for the scalars (see, e.g., [9]). When {2eAQ) = mn, 
the scalar field condenses. We look at a solution with no boundary charge, for which 
{Ao) = 0, i^s = ruH and 

Hence, we're looking at the effects that are l/mn suppressed. 

3. Static Potential Without Fermion Fluctuations 



The uniform background solution obtained above sets a preferred Lorentz frame. 
Let us look at small perturbations in this frame. We introduce perturbations of the 
gauge. An = bf^, and the scalar field, as follows: 

a{x) — ac + t{x) . (15) 

Let us for the time being ignore fermion fluctuations. Physically this would cor- 
respond to the case when fermions are "frozen in" (for instance if they form a 
rigid-enough crystalline structure). We will relax this conditions and include their 
fluctuation in the next Section. 

Then, the Lagrangian density for the fluctuations in the quadratic approximation 
reads [1] 

A = -^41 + lid^rf + ^m'X + 2m^m^ boT . (16) 



Here 



™;.(2e,-^^. (17) 

fniy denotes the field strength for fe^, and we left out for simplicity all the fermionic 
terms. It is useful to integrate out the r field from (16). The remaining Lagrangian 
takes the form: 

A = -i/^„ + iK*; + i;.o^^^*.. (18) 

This Lagrangian contains four components of b^, and no other fields. The first two 
terms in (18) are those of a usual massive photon with three degrees of freedom. The 
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last term is unusual, as it gives rise to the dynamics to the timelike component of 
the gauge field. This term emerged due to the mixing of bo with the dynamical field 
T in (16), and since we integrated out r, 60 inherited its dynamics in a seemingly 
nonlocal way. The fact that there are no pathologies in (16) or (18), such as ghost 
and/or tachyons, can be seen by calculating the Hamiltonian density from (16): 



^^T + T + 2^ + 2 ■ 

Here, ttj = —foj and Pr = Oqt, and the Hamiltonian is positive semi-definite. 



(19) 



The form of the Lagrangian (18) is most useful for calculating a propagator. 
Indeed, the inverse of the quadratic operator that appears in (18) has poles which 
describe all the four propagating degrees of freedom. This propagator, sandwiched 
between two conserved currents and J^, takes the form: 



-P 



mf 



4M^ 

_p2 



(20) 



This describes two transverse photons with the mass m-y, one heavy mode with mass 
2m H, and a light phonon (for their dispersion relations, without and with fermion 
fiuctuations, see Rcfs. [1, 2] respectively). Here a = 1 — AM^uo"^ / p'^m^,{jp' — m^). 
In particular, we are interested in a static potential, which takes the form 



G{u = 0,k) 



4M^ 
1^ 



(21) 



where, we have used the notation M = ^niHiTij, and the first, second, and third 
terms on the r.h.s. of (21) are due to the respective terms in (18). 

Interestingly, when mn ^ Tn^, as is the case here, there is no scale at which 
the photon mass term in (21) would dominate: for > the mass term is 
sub-dominant to the k^ term, while for k^ ^ it is sub-dominant to the Af/k^ 
term. 



Let us rewrite the Green's function (21) in the following form: 



G'(k,u; = 0) = (k2 + m^)"'- 



k^(k^ + m^)^ 
4M4 



(22) 



The first term on the r.h.s. of (22) could be thought of as an instantaneous repul- 
sive screened Coulomb (Yukawa) potential, while the second term as an attractive 
potential due to the phonon. The phonon here is somewhat peculiar - it mixes with 
the timelike component of the gauge field, and thus acquires an instantaneous piece, 
a part of which cancels the instantaneous gauge potential. This is reflected in (22) 
as exact cancellation of the residues of the pole at k^ -|- rriQ — 0. The remaining 
instantaneous part is what's given in (21). 
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As we have already mentioned in Section 1, the remaining two poles of the 
second term in (22), describe respectively the heavy state of mass 2mH, which is 
unimportant for the low-energy dynamics, and a light state - the phonon [2] . Just to 
reiterate, for simplicity of discussions, in this work we're using a somewhat imprecise 
language and calling the whole second term in (22) the phonon contribution. 

The key observation is that at scales larger than 1/M, which are of the primary 
interest, the phonon potential cancels the gauge potential with a high accuracy. Note 
that this cancellation is rehable at scales that are much greater than J^^^^ » M~^, 
and takes place already at scales that are much shorter that the photon Compton 
wavelength 3> Jq 

In a Lorentz-invariant theory having a negative sign in front of a propagator, 
such as the one in the second term in (22), would suggest the presence of a ghost- 
like state. However, this is not the case in a Lorentz-violating theory described by 
our Lagrangian (16) or (18). As we pointed out above, the Hamiltonian density of 
this theory is positive semi-definite, and hence, no ghost or tachyons are present. 
Moreover, consistently with ones expectation, the second term in (22) disappears in 
the limit M — > 0, where Lorentz invariance of (16) or (18) is restored. 

The above described properties can also be seen by calculating the coordinate 
space potential (see, [2]): 

V = {Q.eQ^e) I ^(^^-G{u = 0, k) oc ^i^H^^^^i^ cos (Mr) , (23) 

in which we assumed that r = |x| ^ 1/M and <^ M. Such a potential was first 
derived in [2]. It is sign- indefinite and undergoes modulated oscillations between 
repulsion and attraction. There are an infinite number of points in the position 
space where the force between classical charges would vanish. These are points 
where 

dV 

— (r = rj = 0, n = l,2,... (24) 

The potential wells described by (24) are too shallow to introduce a significant 
quantum-dynamics for static impurities - any two charged probe particles (e.g., 
nuclear hydrogen or helium-3 impurities) separated by a distance r„ ^ Jg , where 
our calculations are reliable, would stay in a static equilibrium as long as K(r„) < 0. 
This is because their charges will be screened efficiently by the ambient condensate. 

Before turning to the next section we make three important comments: 
Comment 1: Concerns the decoupling of the helium-4 sates in the niH oo 
limit. As was already mentioned in Sections 1 and 2, the decoupling takes place 
via "freezing" out of the phonon associated with the collective motion of the heavy 
modes. The fact that this indeed is the case can be seen by looking at the dispersion 
relation for the phonon which is obtained from the poles of the Green's function (21) 
in which nonzero uu is restored. In Refs. [1] and [2] the full dispersion relations for the 
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phonon without and with fermion fluctuations was given respectively. The "freezing 
out" of the phonon takes place in both cases. 

For simplicity we discuss here the case of [2] : For the relevant momenta ^ 
the dispersion relation reduces to a; ~ m-y (1 + k2(k2 - m2)/8M^), and the phonon 
group velocity 

m^|k|(2k2-m2) 

4M^ ' ' (25) 

vanishes in the ttih oo limit. 

Note that for k^ ~ m^/2 the phonon group velocity vanishes for finite mn- This 
describes a state of a nonzero momentum but zero group velocity. The energy of this 
state is also nonzero, and to a good approximation equals to m^. These properties 
are similar to those of a roton in superfluid helium II. Moreover, for excitations with 
k^ > m^/2 the group velocity is positive, while in the opposite case, k^ < m^/2, it 
becomes negative (i.e., the direction of the momentum and that of group velocity 
are opposite to each other). These excitations resemble the positive and negative 
group velocity rotons in superfluid helium II. 

Comment 2: Concerns the issue of applicability of the linearized approximation 
(16) that we adopted in this work. For this we should restore back all the non-linear 
terms that were ignored in (16) and compare them with the linear ones retained 
there. Although this can be done in full generahty, it is instructive to look at the 
dynamics of a longitudinal polarization of the massive vector fleld (i.e., the would-be 
Nambu-Goldstone boson) which determines the fastest growing non-linear terms in 
the Lagrangian at high energies [16, 17]. This is easily achieved by employing the 
substitution — > {d^n)/m^, and taking the decoupling limit: — > and e — > 0, 
while keeping the ratio v = m-y/e flnite and flxed. The resulting Lagrangian for the 
longitudinal mode n, and the remaining fleld r takes the form 

(g^ + + + 2m,(don)r + "^"(^°"^"' + . (26) 

2 V 2v^ V 2 

From this we deduce that the non-linear terms become comparable with the linear 
ones when r ~ we also require that fluctuations of gauge fleld to be smaller 
tat the scalar chemical potential, which gives don ~ mnv. Thus, the domain of 
applicability of the linearized results of this work is 

/J \V2 

T<^;=f^j , don <^ niHV ^ {Jomnf^^ ■ (27) 

In the limit rriH — ^ oo, the domain of applicability of the linearized results shrinks 
to zero. This suggests that the geometric size of the region in which one can mean- 
ingfully talk about the charged condensate should be greater than a certain critical 
size that scales as ~ v~^. The latter tends to inflnity as mn — > oo. 
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Comment 3: The dynamics of mall fluctuations can also be obtained from the 
effective Lagrangian (8). For this we restrict ourselves to the quadratic order in the 
expansion of V{x) 

= 2eE^5o - (^^■^)^ + (2^)^^^^^ + C.(2e)^E^5o^ , (28) 

ZTTlff 

and expand the above expression (amended with the gauge-kinetic and fermionic 
terms) to the second order in the following fluctuations: 

S = \I^+V^t{x) , B^ = h^. (29) 
The resulting Lagrangian for the fluctuations reads: 

A = -\flv - \i^3^f + ^^0^0 - ^^^i + Sm^m^ 6oT , (30) 

where rn^ = 2e^C2Jo ^^^^ before, = (2e)^ Jo/2mj7. This should be compared 
with (16). We notice that unlike in (16), there is no kinetic term for r in (30). This 
is because we decoupled the time-dependence of the heavy field in (6) in the non- 
relativistic limit; thus, the E field, or its fluctuation r, are left to be instantaneous. 
This, obviously, does not change the conclusions on the static potential discussed in 
this Section. 

Furthermore, depending on the coefficient C2 the "electric" and "magnetic" 
masses of the gauge field can be different. The value of C2 cannot be fixed within the 
effective Lagrangian approach. One can only impose a bound on it. For instance, 
without the bo — t mixing, superluminal group velocity is avoided when ^ 

> . (31) 

In this Section the value of C2 was assumed to saturate the above inequality, as 
we put rriQ = m^. However, the fermion fluctuations, which we ignored so far, will 
produce a large hierarchy between the coefficients of the 6g and 6| terms in the 
Lagrangian, even if these were assumed to be equal in the classical theory. 

Finally, let us look at the decoupling limit: e — > 0, while mo/e and m^/e are 
finite and fixed. Then, the Lagrangian (30) reduces to 

777^ 1 1 

^{donf - -^{d,nf + 2mH{don)T - -{d.rf . (32) 



^Notc that the longitudinal part of hj is a propagating mode. Imposing transversality on hj 
would contradict general equations of motion. Also, superluminal group velocity in the present case 
would suggest some instability of the background. In general, superluminal group velocities can 
exist in highly absorptive media, in which case group velocity does not approximate the velocity 
of propagation of a signal. 
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This should be compared with the quadratic part of (32). Again, there is missing 
kinetic term for r, for reasons explained above. As before, r can be integrated out. 
The resulting dispersion relation for the remaining mode reads 

2 

" Am]j + k2(m2/m2) ' ^^^^ 

which for small momenta, <C m^u{m^ / 171^) , gives uj ~ Is? /2mH-i as it should be 
the case for a collection of free massive particles. 



4. Quantum Screening of Static Charges 



In the previous section the fermions were treated as "frozen". In may physical 
circumstances, and in particular in white dwarfs, this is not a good approximation. 
For the helium-electron system the fermion fluctuations should be taken into ac- 
count. This was done in Ref. [2] using the Thomas- Fermi (TF) approximation. The 
result is that the timelike component of the gauge field acquired an additional mass 
- the electric mass. This breaks the degeneracy between the coefficients of 6q and 
6| discussed in the previous Section. 

However, the TF approximation does not capture one significant property of the 
fermion system related to a possibility of exciting gap-less modes near the Fermi 
surface, which we will incorporate below into our results by using the one-loop 
correction to the propagator of (18). For this, we restore back in the Lagrangian 
(18) the fermion kinetic, mass and chemical potential terms and, upon calculating 
the gauge boson propagator, will take into account the known one-loop gauge boson 
polarization diagram. 

Since we're interested in a static potential, we look at the {00} component of 
the propagator Dqo, and the static potential obtained from it: 

V{k) = -Dooiu; = 0, k) = ( k^ + + ^ + F(A;^ fc^, m/) j . (34) 

where the function F^k"^, kF,mf ), which is due to the one-loop photon polarization 
diagram, includes both the vacuum and fermion matter contributions {kp denotes 
the Fermi momentum). A complete expression for Fik"^, kp-iTLf) can be found in Ref. 
[5]. We concentrate on the expression for F{k'^,kF,mf) in the massless (m/ = 0) 
limit that is a good approximation for ultra-relativistic fermions: 

Here /^o stands for the normalization point that appears in the one-loop vacuum 
polarization diagram calculation. The function F introduces a shift of the pole in 
the propagator, corresponding to the "electric mass" of the photon. This part of 
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the pole can be incorporated via the TF approximation, as it was done in [2]. In 
addition, however, the function F also gives rise to branch cuts in the complex |k| 
plane (see [3] for the list of earlier references on this). 

As in the previous section we can decompose the static potential as follows: 

T/(k, ^ 0) ^ (k^ + + F)-' - (^k^ + + F + + ^^' ^ \ (36) 

Note that the first term in (36) is just the instantaneous screened-Coulomb (Yukawa) 
potential of a massive photon with the one-loop polarization correction. Our main 
interest is at distances smaller than m^^. A sphere of radius m^^^ encloses many 

particles within its volume since ^ Jg At these scales, the first term in 
(36) can be approximated by: 

The above expression has a regular pole corresponding to the acquired "electric" 
mass of the photon due to the polarization diagram. The contribution of this pole 
would give rise to an exponentially decaying potential e~"*'^'''/r, where rriei ~ e///. 
This is just an ordinary Debye screening. 

However, as was mentioned above, the expression (37) also has branch cuts in the 
complex |k| plane for k = -^2kp. These branch cuts give rise to the additional terms 
in the static potential which are not exponentially suppressed, but instead have 
an oscillatory behavior with a power-like decaying envelope. In a non-relativistic 
theory they're known as the Priedel oscillations [3]. In the relativistic theory they 
were calculated in Refs. [4, 5] (wc follow here [5] and for simplicity ignore the 
running of the coupling constant due to the vacuum loop) : 



Qigsttem sin(2A;Fr) 



= ,\ . ■ (38) 

47r klr"^ ^ ^ 



Note that these branch cuts have a physical interpretation: Since there is no mass 
gap in the fermion spectrum, a photon can produce a near-the- Fermi-surface particle- 
hole pair of an arbitrarily small energy and the momentum close to ±2^^;'. The imag- 
inary part of the one-loop photon polarization diagram should include the continuum 
of such near-the-Fermi-surface pairs. These are refiected as logarithmic branch cuts 
in the expression for F. 

Thus, if the phonon term (the second term) on the r.h.s. of (36) were absent 
one would have a power-like behavior (38) of the static potential at scales r ^ m~^. 
The phonon term, however, significantly reduces the strength of this potential. The 
result for it could be calculated by directly taking Fourier transform of (34). The 
dominant contribution comes from the branch cuts at A; = ±2/c^. Drawing the 
contours around these cuts in the upper half plane of complex |k| [3, 5], one deduces 
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the result. In the approximation M ^ kp ^ rn^, which is relevant for our system, 
a static potential between like charges scales as 

This potential has an additional suppression factor of 16{kF/M)^ ~ 10~^'^, as com- 
pared to AV in a theory without phonons (38). Nevertheless, AV in (39) dominates 
over the exponentially suppressed part of the total potential found in Ref. [2] and 
discussed in the previous section (see, (23)), for separations between probe particles 
large-enough for the effective field theory description to be applicable. The net static 
interaction in the charged condensate, set by (39), is very weak. It is, however, still 
much stronger than gravitational interaction between a pair of light nuclei. 

Although formally AV in (39) is proportional to a^^, to a good approximation 
it is independent of ctem since oc e^(mj^ Jq)^. 

As in the previous section, the potential for the static charges (39) is shallow 
and sign-indefinite. There are, however, an infinite number of position-points where 
the force between static charges would vanish. These are determined by extremizing 
(39) and obtaining 

4tan(2A;Fr„) = 2A;j7r„, n = 1, 2, 3, ... (40) 

Some of the r„'s are local minima, and thus, static charges (e.g., helium-4 nuclei, 
nuclear hydrogen or helium-3 impurities) placed in those minima with r„ > Jq 
where our approximations are valid, would stay in a local equilibrium due to efficient 
screening by the condensate. Hence, the condensate could tolerate some fraction of 
charged impurities before the latter could significantly affect the condensate itself. 

Finite temperature corrections would modify these results quantitatively, how- 
ever, the main point of efficient screening of the static charges in the condensate 
should remain valid at temperatures well-below the condensation point. For in- 
stance, in white dwarfs with temperature 10^ — 10^ K we would expect the dom- 
inant temperature-dependent corrections to the potential to be proportional to 
T/Jq^ - (10-^ - 10-3) < 1, which are negligible. 
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